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We explicitly relate effective meson-baryon Lagrangian models, chiral bags, and Skyrmions 
in the following way. First, effective Lagrangians are constructed in a manner consistent 
with an underlying large- A'^c QCD. An infinite set of graphs dress the bare Yukawa couplings 
at leading order in 1/Nc, and are summed using semiclassical techniques. What emerges is 
a picture of the large- A^c baryon reminiscent of the chiral bag: hedgehog pions for r > A~^ 
patched onto bare nucleon degrees of freedom for r < A~^, where the "bag radius" A~^ 
is the UV cutoff on the graphs. Next, a novel renormahzation group (RG) is derived, in 
which the bare Yukawa couplings, baryon masses and hyperfine baryon mass splittings run 
with A. Finally, this RG flow is shown to act as a filter on the renormalized Lagrangian 
parameters: when they are fine-tuned to obey Skyrme-model relations the continuum limit 
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1. Introduction 

1.1. Effective hadron Lagrangians versus Skyrmions 

In the absence of reliable quantitative methods for computing the low-energy proper- 
ties of QCD, a wide variety of phenomenological models of the nucleon have emerged and 
flourished. One general class of models, which predates QCD by some 30 years, starts from 
an effective Lagrangian for the baryon and meson fields. The hope in this approach is that 
the relevant physics is contained in the complete set of hadron Feynman diagrams. An 
orthogonal approach, pioneered by Skyrme in the early 1960s, uses topology: the baryon 
is viewed as a soliton, or Skyrmion, in the field of mesons [|l],|2|. 

On the face of it, the two approaches could not be more opposite. Baryon number in 
an effective quantum field theory of hadrons is simply the Noether charge associated with 
a [/(I) symmetry of the Lagrangian. In contrast, in the Skyrmion picture, baryon number 
is not associated with any continuous symmetry but is instead a topological invariant: the 
winding number of a meson field configuration. And the physics of Skyrmions is expressed, 
not in the language of Feynman diagrams, but in vocabulary more appropriate to solitons 
and other extended objects: collective coordinate quantization, symmetry classification of 
small fiuctuations about the soliton, and so forth. 

In this paper, we exhibit a precise connection between these two disparate approaches 
(see Fig. 1). The bridge between them is built by combining in a new way two important 
theoretical constructs: 't Hooft's 1/Nc expansion Nc being the number of colors in the 
underlying gauge theory, and Wilson's renormalization group 0. The parameter governing 
the fiow of this so-called large-Nc renormalization group is an ultraviolet cutoff A ~ A^^ 
which regulates the divergences in the effective Feynman diagrams. Our main result is that 
soliton models and (suitably fine-tuned) effective Lagrangian models are completely equiv- 
alent at leading order ,0 in the "continuum limit" A ^ oo (to borrow lattice terminology). 
In this Section we introduce the key ideas behind the large- A^ renormalization group; the 
remainder of the paper is primarily devoted to exhibiting its solutions in a series of simple 
models. 

^ Prom now on, "leading order" refers to the 1/Nc expansion. The need to fine-tune the 
renormalized Lagrangian couplings is explained in Sec. 1.4 below. 
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1.2. The large-Nc limit and semiclassical physics 

In 't Hooft's original formulation, the 1/Nc expansion is studied directly in QCD; as 
Nc ^ oo the physics is dominated by the quenched planar quark-gluon graphs. But the 
large- A^c limit is actually much more predictive when the lessons of planar QCD are imple- 
mented, not at the fundamental quark-gluon level, but rather at the level of phenomeno- 
logical models such as above. Thus, on the one hand, in effective hadron Lagrangians, it 
implies a set of large-Nc selection rules which place powerful constraints on the allowed 
particle spectrum and couplings contained in one's model. As examples of such rules 
whereas meson masses typically scale like A^^, baryons (being made up of Nc quarks) have 
masses that grow linearly with Nc. And whereas Yukawa couplings are strong, growing like 



\/iVc, meson self-interactions are weak, with n-meson vertices disappearing like Nc 

The complete set of these large- A'^^ selection rules is reviewed at the beginning of Sec. 2 
below. On the other hand, in Skyrmion physics, 1/Nc plays an ostensibly different role, 
which is more easily summarized: it parametrizes the semiclassical expansion about the 
soliton. This is because 1/Nc always enters into Skyrmion Lagrangians in the combination 



In fact, this difference is illusory: the effective Lagrangian approach, too, becomes 
semiclassical in the large- A'^c limit (a key to the equivalence of the two pictures). What 
we mean by this is twofold. First, in calculating the leading-order contribution to meson- 
baryon Green's functions, the naive, graph-by-graph, perturbative method fails, and one is 
forced instead to sum an infinite class of diagrams. Second, this sum may be accomplished 
by solving classical equations of motion for the meson fields in the background of the 
baryon source. 

In order that these two points be understood, let us be very explicit at this stage, 
and in so doing, introduce the central field-theoretic problem of this paper. Consider the 
bare Yukawa coupling (/bare depicted in Fig. 2a, which scales like y/Nc as stated above. In 
principle, we would like to sum all radiative corrections to this vertex such as Figs. 2b-d, 
and thereby extract the renormalized Yukawa coupling Qren, to leading order in l/A^c-^ 
Focus first on Fig. 2b. It contains two factors of (/bare, therefore two factors of ^/N^, as 
well as a 3-meson vertex which goes like 1 / ^/N^; so this graph too scales like ^/N^, and is a 

^ The calculation of renormalized single- meson emission/absorption from the nucleon is the 
simplest arena for our semiclassical methods, which may also be applied to more complicated 
processes such as meson-baryon and baryon-baryon scattering 




h/Nc. 
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leading- order correction to (/bare- Likewise Fig. 2c scales like ^/Nc, as the reader can check 
by multiplying all vertices together. An example of a subleading correction is Fig. 2d. 
Unlike the others it contains a purely mesonic loop (indicated by the arrow), hence two 
extra factors of 1/ y/N^ uncompensated by extra Yukawa vertices. So Fig. 2d and the like 
are 1/Nc corrections, and will not concern us further. 

A moment's thought confirms the rule: the leading-order dressings are the infinite set 
of diagrams which contain only meson trees if the baryon line is erased. And the sum of 
all trees — like a soliton — is given by the solution to a classical equation of motion, which 
we write down in Sec. 3 below. The role of the bare large- A''^ baryon in this equation is 
that of a heavy, slow-moving source, smeared out over a length-scale A~^ so as to cut 
off the short-distance divergences in the original Feynman graphs. As we shall see, the 
solution of the regulated equation is a hedgehog cloud of pions and/or other mesons for 
r > A~^, glued onto the bare nucleon degrees of freedom which are restricted to r < A~^. 
The energy of this cloud renormalizes the mass of the baryon (Mbare -^ren) while its 
large-distance behavior determines the physical Yukawa coupling ((/bare dren)- 

1.3. Of chiral bags and Cheshire cats 

The resulting picture of the meson-dressed large- A'^c baryon is highly reminiscent of 
yet a third class of phenomenological models, the chiral bag models [p!0H T3|. These, too, are 
hybrid descriptions of the dressed baryon, in which explicit quark (rather than nucleon) 
degrees of freedom inside a bag of radius R are matched onto an effective theory of hedgehog 
pions outside the bag. Even this presumably important distinction between 'nucleon' 
versus 'quark' degrees of freedom inside the bag disappears as Nc ^ oo.i For, in this 
limit, the N^. quarks may be treated in Hartree approximation, and their individual wave 
functions effectively condense into a common mean-field wave function, which we may 
identify with the "wave function of the nucleon." Outside the bag, the analogy is closer 
still: the pion field configuration is again determined by solving a nonlinear field equation 
coupled to a static source a,t r = R. The only significant difference between our composite 
meson-dressed large- A^ baryon and the traditional chiral bag is this: our composite baryon 

^ This observation was originally made by Witten (Ref. Q, Sees. 5 and 9), and exploited by 
Gervais and Sakita (Ref. [Q, Sec. V). These two papers are highly recommended background 
reading, as they too are concerned primarily with the semiclassical nature of large- A^c- In partic- 
ular Gervais and Sakita were the first to study chiral-bag-type structures in this limit, although 
not from our effective hadron Lagrangian starting point. 
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follows solely from large Nc and has nothing whatsoever to do with chiral symmetry! For 
this reason we shall refer to it as a "chiral bag" (Fig. 1 again), being careful to retain the 
quotation marks to avoid confusion with the traditional chiral bag. 

Associated with the traditional chiral bag in the recent literature is the so-called 
Cheshire cat principle [|T5| , p!3[| , which states that physical quantities should be independent 
of the size and shape of the bag. This proposal is motivated by the success of bosonization 
in 1+1 dimensions, which provides an exact mapping between, for instance, the elementary 
fermion of the massive Thirring model and the soliton of the sine-Gordon equation ||16|| . 
In its most extreme form, it implies that R can safely be set to zero, yielding a description 
of the nucleon purely in terms of pion fields. One may conjecture that this limit is nothing 



other than the Skyrme model [p!2| , p!3| . However, since no equivalent of bosonization has 
been made to work in 3+1 dimensions, the theoretical status of the Cheshire cat principle 
has remained uncertain. 

From our present large-A''^ perspective, the bag radius R is best reinterpreted as 
the short-distance cutoff of the large- A'^c effective theory. It is then easy to see the 
significance of the Cheshire cat principle: the cutoff-independence of physical quantities 
is usually referred to as renormalization group invariance Q]. However, except in special 
cases, renormalization group invariance of the physical masses and couplings never comes 
for free; the corresponding bare quantities must be varied simultaneously with the cutoff. 
As stated earlier, we will refer to this program as the large-Nc renormalization group, 
and devote considerable space to mapping out its solutions. To our understanding, the 
concept of such a flow is not only absent from the usual Cheshire cat philosophy, but in 
fact orthogonal to it; in that scheme is nothing more than a gauge-fixing parameter. 



1.4- The large-Nc renormalization group as a "filter" 

The limit of zero bag size, ^ or A ^ oo, corresponds to removing the cutoff 
completely. This continuum limit, if it can be taken at all, is by definition a UV stable 
fixed point of the renormalization group fiow.@ In light of the above discussion, it is 
natural to conjecture that such a fixed point exists if and only if the homogeneous meson 

^ Because we are implementing renormalization group invariance strictly at leading order in 
1/Nc, the limit A — > oo is to be taken after Nc — > oo. We suspect that these limits do not 
commute. 
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field equations {i.e., with the baryonic source set to zero) supports a sohton solution.El 
Thus we propose that variants of the Skyrme model describe the possible continuum limits 
of effective Lagrangian theories of mesons and baryons in the large- A''^ limit. 

On reflection, there is an obvious counting problem with this scenario: effective La- 
grangians always contain more free parameters than the corresponding soliton models. 
Thus, in the former, the physical masses and couplings are all independent, while in the 
latter, there exist non-trivial relations among them; for instance, the Yukawa constant 
Qren is Completely determined by the meson self-couplings [|2|. (This feature, of course, is 
precisely the point of the Skyrmion approach.) If our proposal is correct, it follows that 
unless the renormalized parameters in the effective Lagrangian are tuned exactly to those of 
the corresponding soliton model, there must be some mathematical obstacle to taking the 
continuum limit. In other words, we conjecture that the large- A''^ renormalization group 
acts as a filter, blocking the path to the continuum except for a measure-zero subset of 
the space of renormalized parameters. This fllter idea is the central theme of this paper, 
and is explicitly realized in the models to follow. The attentive reader will recognize in the 
phrasing of this counting problem the classic symptom of the existence of an "irrelevant 
operator" [|| ; naturally this operator turns out to be the bare Yukawa coupling itself. 

1.5. The plan of this paper 

The paper is organized as follows. Sections 2 and 3 recapitulate our recent Letter 
in somewhat more detail, and at a more leisurely pace. (Independently, Manohar has 
reached similar conclusions to Ref. 0; other relevant precursors are Refs. and [|n| .) 
In particular. Sec. 2 reviews the large- A'^c selection rules mentioned above, and gives the 
recipe for constructing large- A^c-compatible effective Lagrangians. Section 3 is devoted to 
semiclassics: the problem posed in Fig. 2 is solved completely at a formal level, the issue 
of meson contributions to the baryon self-energy is examined, the hedgehog structure of 
the meson cloud is revealed, and the large- A'^c renormalization group is deflned. 

The reader already familiar with the contents of Ref. is encouraged to skip di- 
rectly to Sees. 4-8, in which the large-A'c renormalization group is applied to a series 
of effective hadron models. What distinguishes these models from one another is our 
choice of the purely mesonic piece of the action, £meson- In Sec. 4, £meson is simply the 

^ The only exception we have found to this rule is a theory of non-self-interacting mesons 
(Sec. 4). 
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free pion Lagrangian. Not surprisingly, one finds essentially no running of the bare nN 
Yukawa coupling: (7bare(A) = (/ren- Much more interesting is the model of Sec. 5, in which 
C-meson = ^ Tt d^U^d^U , the leading term in chiral perturbation theory. This is a clean 
initial test of the "filter" idea conjectured above, because this model is known not to sup- 



port a soliton, thanks to Derrick's theorem |]T^. And indeed, we discover a critical value 
of the cutoff, Acrit — 340 MeV, beyond which the large- A^c renormalization group cannot 
be pushed. In Sec. 6 we construct an analytically soluble 2+1 dimensional model, in which 
the would-be Skyrmion is simply the instanton of the 0(3) a model in one lower dimen- 
sion. In this model the physics is all that one might have hoped: unlike Sec. 5 there is no 
ultraviolet obstruction to the large- A'^c renormalization group; the bare coupling (/bare (A) 
is irrelevant, flowing to zero like A~^; and a "toy" Skyrme model indeed lies at the end 
of the large- A'^c renormalization group trajectory, at which point explicit baryon number 
effectively transmutes into the winding number of the meson cloud. 

Finally, in Sees. 7-8 we augment the nonlinear pion Lagrangian of Sec. 5 with the 
4-derivative "Skyrme term" and study the large- A'^c renormalization group, both in the 
measure-zero case when the physical Yukawa coupling is tuned precisely to its Skyrme- 
model value (Sec. 7), and in the generic case when they differ (Sec. 8). In a surprising way, 
involving essential singularities in Skyrme's equation, and local instabilities that develop in 
the pion cloud for A > Acrit (see the Appendix for technical details), our "fllter" conjecture 
is borne out. 

In light of this already long Introduction, we will spare the reader a Conclusions 
section, and commend him instead to keep the "large- A'c renormalization group as filter" 
idea firmly in mind as he works his way through the examples. 



2. Constructing large-A'c effective Lagrangians 

2.1. Large-Nc selection rules 

As mentioned in Sec. 1, the large- A^c limit imposes several stringent requirements 
on the allowed spectrum and interactions of hadrons, which we now review. These rules 
can be derived quite independently from several different approaches: not just a planar- 
diagrammatic analysis of large-A^c QCD [[7|,|8| JT9| , pO| but also the Hartree approximation 
originally employed by Witten 0, the Skyrme model [p|, p^ , |2T[] , the non-relativistic quark 
model [p2| -p4|, and finally the self-consistency of the effective hadron Lagrangian [p^p5| 
as recently emphasized in a series of interesting papers by Dashen, Jenkins and Manohar 
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|2^,27|. The following rules (the first two of which were already invoked in Sec. 1) should 
be considered robust, model-independent features of the large- A^c limit: 

(i) As noted originally by Veneziano, purely mesonic vertices with n external legs 
scale like Nc An important example with n = 1 is that the pion decay constant 
/tt ~ whereas meson masses (n = 2) generically scale like A^^. 

(ii) Vertices with two baryon legs and n meson legs also scale like Nc "'^^ so that 
baryon masses {n = 0) and Yukawa couplings {n = 1) grow like Nc and v^iV^, respectively 

(iii) The 2-fiavor stable baryon spectrum of large- A''^ QCD (with Nc odd) consists of 
an infinite tower of positive parity states with I = J = |, |, To leading order these 
states are degenerate [p|, p^ , |T9| , |20| , p6l , with bare mass Mb 

are ~ -^c ('^•S baryons are made of 
Nc quarks). Hyperfine baryon mass splittings have the form J( J-(- l)/2Xbare where Xbare ~ 
Nc P, |l9ip0| , |27[] . In the Skyrme model these states correspond to (iso)rotational excitations 
of the static hedgehog and the splitting term corresponds exactly to the rotational kinetic 
energy of the Skyrmion [Q. 

(iv) In contrast, the spatial extent of a baryon does not grow but has a smooth Nc- 
independent limit as Nc oo as do the various baryon form factors such as electromagnetic 
charge distributions, and baryon number density itself 

(v) Yukawa couplings are constrained to obey the "proportionality rule" which fixes 
the interaction strength of a given meson with each member of the baryon tower as a 
multiple of one overall coupling constant P, p^pl| , [23| , p6| (e.g., g^NN oc g^NA oc g^AA ■ ■ ■ )' 
up to corrections of order l/N^ P^ . 

(vi) Finally, the allowed couplings of mesons to the baryon tower must obey the 
It = Jt rule pT| , p3| , |25[| . For example the p meson must be tensor-coupled to the nucleon 
while the u meson is vector-coupled at leading order in [p5| , |20| , p4[ , in good agreement 
with phenomenology . When crossed from the t-channel to the s-channel, this rule also 
implies nontrivial model-independent relations among meson-baryon 5" matrix elements 
II29H31I], which are tested against the experimental data in Refs. [|30|-|3l|.B 



^ In the Skyrme-model representation this tower is truly infinite whereas in the quark-model 
representation it tops off at Nc/2; see Ref. [24| for a detailed discussion of how to translate 
Skyrme-model operators into quark-model operators and vice versa in light of this difference. 

^ Rules (v) and (vi) are elementary examples of "large- A^'c group theory," the state-of-the-art 
phenomenological predictions of which are summarized in Ref. [p^ . The existence of such group- 
theoretic relations may be traced to the fact that SU{2Nf) spinxflavor symmetry becomes exact 
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These selection rules must be implemented in any large- A'^c-compatible effective hadron 
Lagrangian. What we stressed in our Letter [|| and will review below is the consistency 
of these rules, meaning that once they are incorporated into the bare Lagrangian they 
continue to hold for physical, renormalized quantities as well. Thus, while the baryon 
mass spectrum is renormalized at leading order by the interactions with the mesons, the 
J( J + 1) structure discussed in (iii) is preserved, and so, too, the form of the 0{\/Wc) 
Yukawa couplings of the mesons to the baryon tower dictated by (v) and (vi) . And while 
the "bare" nucleon size (as given by the naive ultraviolet cutoff ~ A*"^) is effectively 
enlarged by the meson cloud, it remains of order A^^ as per (iv). This is because the 
spatial extent of the cloud is dictated by the parameters (i) of the meson Lagrangian, 
for instance or the product of Skyrme- model parameters (eg/jr)"^ as we shall see in 

Sees. 7-8 below. 

In general, an effective meson-baryon Lagrangian is a sum of four parts: 

-^efF -^baryon "l~ -^meson ~l" -^yukawa ~l" -^seagull ■ (2-1) 

In light of the above selection rules, let us discuss, in turn, the proper construction of each 
of these parts. 

2.2. Constructing Lbaryon- 

We start with a relativistic baryon Lagrangian for the tower described in (iii) above: 

N{iYd^, - Mn)N + (higher I = J baryons) (2.2) 

where A^ means (^). We will need to recast this unwieldy infinite sum in a more useful 
form. Since large-A^c baryons are heavy, it is natural to split their propagators in the 
usual way into forwards-in-time (t/-spinor) plus backwards-in-time (l^-spinor) pieces. The 
latter account for the so-called Z-graph contributions to Feynman diagrams, which turn 
out to be subleading as we review momentarily (Sec. 2.5). The remaining time-ordered 
diagrammatics is one in which Z-graphs have been eliminated, and with them the higher 
components of the baryon-antibaryon Fock space. In this way, baryon quantum field theory 
collapses to baryon quantum mechanics. But we can simplify the physics even further. For 



as Nc oo 1 14, 26]. This is one of the two main attractions of large- Ac physics, the other being 



its amenability to a semiclassical treatment. 
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processes involving one baryon only, interacting with an arbitrary number of mesons, it is 
natural to work in or near the baryon's rest frame, in which case, finally. 



^baryon 



-Mbare + ^MfaareX^ + l^^.^TlA^ A + • • • , (2.3) 



the dots denoting 1/Nc corrections. The first two terms on the right-hand side are the 
usual nonrelativistic approximation to the relativistic mass-energy, X(t) being the baryon's 
position. 

The third term, the (iso)rotational kinetic energy, is perhaps less familiar. It repre- 
sents free motion on the baryon's spin/isospin manifold (just as the previous term denotes 
free spatial translation), with A{t) G SU{2) being the baryon's spin/isospin collective co- 
ordinate. The full meaning of this term is revealed in the beautiful path integral identity 



due to Schulman |]32[| : 



/•A(t2)=A2 r 
/ VA{t) expi dt{- Mbare + Xbarel^i^i) 



~A{t2)=A2 

'A{ti)=Ai 

/ ^ / ^ \ I Z ^Z ' ' z I ' 



(2.4) 



where 



MLe = Mbare + " (2-5) 

^-^bare 



In other words, XbareTrA^^A is convenient shorthand for the free propagation of an infinite 
tower of / = J baryons with the mass spectrum of a rigid rotor — exactly as required by 
rule (iii) given above. We will adopt Skyrme-model nomenclature and refer to Xbare as 
the "bare moment of inertia" of the baryon 0. The brackets 

{A\i=;[) = {2J + l)'/'{-y-^^Di'l^^{A^) (2.6) 

in Eq. (|2.4| ) are just the change-of-basis overlaps between the usual spin-isospin baryon 
representation (the nucleons with / = J = |, the A's with / = J = |, etc., with spin and 
isospin z components Sz and Zz), and baryon states \A) sharp instead in the spin/isospin 
collective coordinate. The \A) basis too, while popularized by the Skyrme model 0, is 
useful more generally in large- A^c physics |lT^ , p2| ,|2l 
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2.3. Constructing L 



As for the meson piece of the action, it is best to leave L^eson completely unspecified 
for the time being, subject only to the scaling rule (i) discussed above. Importantly, if one 
then rescales all meson fields by a mass parameter proportional to \/N^ (such as the pion 
decay constant /^r), a multiplicative factor of Nc/h sits in front of the meson action. So 
a leading-order analysis in the 1/Nc expansion is tantamount to a semiclassical {% 0) 
treatment of the mesonic part of the path integral. We will exploit this feature shortly. 

2.4- Constructing Lyukawa- 

Next, consider Lyukawa- For large A^c, the usual pseudovector coupling of the pion to 
the nucleonji 

glT.d^n-N^^^'^TN , (2.7) 

must be augmented by similar couplings to the entire I = J tower of baryons in a manner 
fixed by the proportionality rule (v). As before, the \A) basis for the baryons permits an 
especially compact representation of this set of couplings, namely 

'^alT. ^^"(^) 5^(x') . (2.8) 

Here the primed space-time coordinate x' = (t',x') is the Poincare transformation of x 
into the center-of-mass frame of the baryon, assumed to be moving with fixed velocity X 
relative to the Lab frame; when X <^ 1 so that Lorentz contractions are irrelevant, one 
simply has 

x' ^ X- X(t) . (2.9) 

The rotation matrix D^^\A) is an operator (hence the 'hat' on the A) on the spin and 
isospin quantum numbers of the single-baryon Hilbert space. It is completely specified by 

^ We have absorbed the normally explicit factor of {2Mn)^'^ into the pseudovector coupling 
constant q^nn, which therefore has dimensions of length. The reason that pseudovector coupling 
is far preferable to pseudoscalar coupling in large- physics is that the latter involves an awkward 
cancellation: 7^ couples the 'large' to the 'small' components of the nucleon spinor, the 'small' 
components being down by 1/Nc, which compensates the fact that the pseudosca/ar constant 

3 /2 

grows like Nc as dictated by the Goldberger-Treiman relation. 
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its matrix elements in the conventional spin-isospin baryon basis: 

JSU{2) 



JSU{2) 

[(2 J + 1)(2 J' + 1)] ^'^ 

1 J J' \ / 1 J J' 



To obtain the first equality we have inserted a complete set of \A) states on which D^^^ [A) 



is sharp; the second equality follows from Eq. ( |2.6| ); and the third, from the textbook 
expression for the integral of three Wigner D-functions. The two resulting 3j symbols 
express conservation of isospin and angular momentum, respectively, while the overall 
square-root coefficient embodies the proportionality rule (v). 

It is easily checked that the terms in Eq. ( |2.7|) involving the spatial derivatives of 
the pion are correctly reproduced by Eqs. ( |2.8|) and (|2.1C1| ), once one specializes to 'in' 
and 'out' nucleons by plugging in J = J' = |. In contrast, the time derivative of the 
pion has been dropped in moving from (|2.7| ) to (|2.8| ). This is because d^ii multiplies the 



Dirac matrix 7^7^ which couples the 'large' components of the baryon's Dirac spinor to 
the 'small' components, the latter being down by f/c ~ l/Nc- As often happens, the l/N^. 
expansion has broken apart a Lorentz-invariant quantity. 

Beyond the nucleons, the coupling ( p.8|) contains useful phenomenological information 
about the higher I = J baryons as well. Sandwiching it between nucleon and A states using 
Eq. ( |2.1C1| ), one calculates a A ^ Nn decay width within a few MeV of its experimental 
value of 120 MeV [^|6|,|33[] . In the same way, one discovers that the higher baryons have 
widths so large (Es^^^ ^ 800 MeV, ^r^i+n ~ 2600 MeV, ^^^l+n ~ 6400 MeV, etc.) 
that they cannot sensibly be regarded as "particles" at all [|5|,p3|. So these higher-spin 
states, the existence of which is often considered a major failing of the large- A^c approach, 
are actually not in conflict with phenomenology.! 



The story of these unwanted large-A^c baryons is much more complicated in models with 



three or more light flavors. The 3-flavor Skyrme model |34], for example, predicts exotic baryons 
not just with high spin but also with low spin, for instance a spin-^ antidecuplet. In quark- 
model language such exotics map onto baryons with Nc quarks plus extra qq pairs [22|. Large- A^c 
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Like the pions, the aUowed nucleon couphngs of the p, uj and/or cr mesons (for example) 
must also be echoed by couplings to the entire I = J tower : 



g^'^d^P.-Na^^fN 3<7,^-^e,,,^p«(x)I^W(A(t'))53(x') (2.11a) 

gtT.^.N^N gtT.u^'ix)5'i^') (2.116) 
^^^-aiViV <7^-Xa:)53(x') , (2.11c) 

again dropping Dirac structures that involve 'small' components. The scalings are 

bare ^ bare ^ bare ^ bare ^ / at ('q -i r)\ 

y-KNN ~ iJpNN ~ iJuJNN ~ yaNN ~ V "''c • l^'-l-^J 

Conspicuously absent from this list of permitted Yukawa couplings are the vector cou- 
pling of the p meson and the tensor coupling of the u meson, g^^p^j^^fP^ ■ N'^^tN and 
g^^%d^u!,yNa'^'^ N . These alternative couplings are forbidden by rule (vi) given earlier, 
meaning that in contrast to ( |2.12D , 



iJpNN iJujNN /jV^ \ 



and we can forget about them. The incorporation of these additional mesons is deferred to 
future work;0 for simplicity the explicit meson models analyzed below will be built from 
pions alone. 



counting implies that each such pair costs a factor of 1/ \/Nc to produce [^,^. And indeed one 
can show that matrix elements of physically relevant operators between "normal" and "exotic" 
baryons are down precisely by where k is the number of such qq pairs; see Ref. [ |24| , 

Sec. XIII for details. 

A natural extension of the ideas in the present paper, and one more in keeping with Wilson's 
original philosophy of the renormalization group is to add more and more mesons to the model 
as the cutoff A is increased. This is consistent with the spectral representation of large- A^c physics, 
which requires, not just an infinite tower of baryons, but apparently also an infinite number of 
mesons in each channel The smearing of the baryon over a distance A^^ can then be thought 
of as being due to the interactions with mesons of mass greater than A, which are omitted from 
the model, whereas all mesons of mass less than A are explicitly kept. 



12 



2.5. Constructing Lseaguii 

Finally we discuss couplings such as Fig. 3a, in which more than one meson interacts 
with the baryon at the same space-time point. For want of a better word we refer to 
these vertices as "seagulls." From QCD one can show that n-meson seagulls generically 
scale like Nc (rule (ii) above). A familiar example from chiral perturbation theory is 
the coupling to the nucleon of the pion's axial current, which is formed from the U field, 
U = exp (2z7r ■ t/Ztt). Taylor expanding U gives coefficients proportional to (7^/"" for the 
n-pion couplings to the nucleon, where n = 1, 3, 5, ■ ■ ■. Since (7^ ~ Nc and /^r ~ y/Nc these 
indeed have the full-strength scaling behavior with Nc. An exception to this scaling rule is 
the set of n-pion seagulls arising from the vector current coupling. These are proportional 
to gyf~^, n = 2, 4, 6, ■ ■ ■ , and since gy = Ihy vector current conservation they all drop 
out as Nc 00. 

Quite aside from any "bare" seagulls that one may choose to include in the bare 
Lagrangian, one must also examine the effective seagulls that arise from the backwards- 
in-time baryons (Z-graphs), once a Feynman diagram is decomposed into a sum of time- 
ordered diagrams. Inevitably, these approximately pointlike induced vertices are termed 
"Z-guUs." As illustrated in Figs. 3b-c, they come from approximating the IZ-spinor propa- 
gator by the constant — z/2Mi.en up to nonlocal 0{Nc'^) corrections, where Mj-en is the phys- 
ical baryon mass. Naively, the strength of the Z-guU in Fig. 3c is then (/bareZ-^ren ~ Nc; 
since this is the same order as a full-strength 2-meson bare seagull, it appears such a ver- 
tex must be kept. However, this naive counting is incorrect^ The reason is easily seen 
by writing out the Feynman rules for the U and V spinors directly. For precisely those 
Yukawa couplings ( |2.8| ) and ( |2.11j ) permitted by the It = Jt rule (vi), there is an extra 
suppression of a ■ p/M ~ 1/Nc at each of the two vertices in Fig. 3b. This is the cost of 
turning a [/-spinor into a V^-spinor or vice versa, when the meson coupling is dominantly 
block-diagonal in the Dirac space. We conclude that in self-consistent large- A^c models, 
Z-gulls are actually suppressed by two powers of 1/Nc compared with bare seagulls, and 
may safely be dropped. 

For simplicity, in the explicit models analyzed below, we will choose to set to zero 
all bare seagulls as well. Instead we will focus on the renormalization of the Yukawa 
interaction and show that, in many cases, it corresponds to an irrelevant operator of the 
large- A'c renormalization group. We conjecture that in these cases the higher-point bare 
seagull couplings are also irrelevant operators, but we will leave this interesting question 
for future work. 

We are indebted to Jim Friar for pointing out an error in an earlier draft of this paper. 
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3. Large-A^c Semiclassical Analysis 

3.1. Formal summation of the leading- order graphs 

We now return to the central problem posed in Sec. 1, namely the summation of 
the leading-order contributions to the renormalized Yukawa constants Qren- Recall from 
Fig. 2 that the leading-order graphs are those containing no purely mesonic loops, in other 
words, those graphs which would be meson trees if one were to erase the baryon line. The 
complete set of such graphs is captured in Fig. 4. 

Not surprisingly, being tree-like. Fig. 4 can be generated as the solution to a classical 
equation of motion. As a concrete example, with pions only, suppose that 

i^meson = = ^ {8 ^n)^ - ^mln^ - V (n) (3.1) 

where V contains the quartic and higher pion self-interactions. The pion trees with one 
external pion line sum to a quantity we call tTci (the subscript 'cl' standing interchangeably 
for 'classical' or 'cloud') which solves the Euler-Lagrange equation implied by Eqs. ( |3.1J ) 
and (13): 

(□ + m^)7rS(x,t)-F^ = 3glT.Dli\A{t'))£^S'{^') , (3.2) 

x' referring to the center-of-mass frame of the moving baryon as before. This equation is 
illustrated in diagrammatic form in Fig. 5. Comparing Fig. 5 with Fig. 4, we find apparent 
agreement, save for the "missing" sum over the n! tanglings. But this sum is already 
implicit in Fig. 5. To see this, insert into the right-hand side of Fig. 5 the resolution of 
unity into a sum over all n\ (boosted) time orderings of the attachment points Zk, 

1 = E ^«(2) - <(1))^«(3) - <(2)) X ■ ■ ■ X H^X) - ' (3-3) 

p being a permutation, and for each element in this sum relabel 2p(fc) ^ ^fc- In this way 
we recapture, not Fig. 4 precisely, but rather Fig. 6, which differs from Fig. 4 only in 
the time-ordering prescription up the baryon line. This difference is truly unimportant; 
backwards-in-time baryon propagation is always a 1/Nc effect .0 



In this particular instance the suppression of Z-graphs is even greater than that discussed 
in Sec. 2.5 above; shrinking a y-spinor propagator to a point here results in a meson loop, which 
itself is 1/Nc suppressed (see Fig. 2d). 
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To summarize, Eq. (|3.2|) correctly accounts for the leading-order contributions to 
the renormalized meson-baryon coupling, up to 1/Nc corrections. The answer is conve- 
niently expressed as a quantum mechanical path integral over the baryon's translational 
and (iso)rotational collective coordinates: 

j 'DA{t)'D^{t) 7r"i(x, t) expi j (it (^Lbaryon + -^meson[7rcl] + i^yukawa[7rcl] j , (3.4) 

ignoring seagulls for simplicity as stated above. Our progress to this point has been that 
the path integration over the field-theoretic variable tt has been carried out in semiclassical 
approximation. This simply means replacing tt by tTci everywhere in ( ^.4] ), the justification 
being that 1/Nc always appears in the combination fi/Nc as highlighted earlier. 

Thanks to the rescaling argument of Sec. 2.3, the meson cloud, like a skyrmion, is a 
nonperturbatively large configuration, scaling like /^r ~ a/ZVc, although its spatial extent 
goes like A*"^. Another general point: it is not sufficient that the cloud be a solution to 
the Euler-Lagrange equation; it must actually be a local minimum, in other words it must 
be locally stable against small deformations. This stability issue will reemerge in Sec. 8 
below. The interesting related question. What does an unstable cloud collapse intdl, is 
examined at the end of the Appendix. 

3.2. Baryon self-energy and baryon-meson vertex corrections 

What is the meaning of the terms i^meson [ttci] +-^yukawa[7rci] in Eq. (|3.4|)? Graphically, 
the answer can be seen in Fig. 7. As Fig. 7c in particular makes clear, these terms 
fully account for the meson-tree baryon self-energy and baryon-meson vertex corrections 
that we have neglected till now [^^. Taylor-expanding the exponential of these terms 
produces an arbitrary number of such insertions at all placements along the baryon line, 
automatically with the correct combinatorics. Moreover, since vfci itself depends on the 
baryon collective coordinates X(t) and A(t) through the Yukawa source on the right-hand 
side of (|3.2|) , Lmeson [ttci] +-^yukawa[7rci] have this dependence as well, and on general grounds 
must have the formlll 

i^meson[vfcl] + Ly^kawa [tTci] = "Md + ^MclX^ + I^iTyA^ A +■■■ . (3.5) 



The third term on the right-hand side, taken together with the identity (2.4), establishes for 
all leading-order baryon self-energy and baryon-meson vertex corrections, and for any value of the 
pion mass, the self-consistency of the hyperfine baryon mass splittings originally noted by Jenkins 



for the simplest such graph in the chiral limit |27]. As for the first two terms on the right-hand 



side, the fact that the same Md appears twice, with a relative weighting of — 4X^, follows trivially 
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Since tTci ~ it follows that Md and Td scale like Nc, just like Mbare and Xbare- The dots 
in Eq. ( |3.5D indicate terms suppressed in the 1/Nc expansion, such as XA cross terms.0 
The lesson of Eq. (|3.5|) is that the classical meson cloud tTci (likewise pd, Wd, C7"d, etc.) 
gives a form-preserving and leading- order^ correction to the 'bare' baryon expressions 
(^)-(|23p, so that effectively 

M^^, — > Mi^ = M,en + ^^^^ (3.6) 

where the renormalized baryon mass and moment of inertia are simply the sumS 

Mren — Affaare + -^cl , ^ren = ^bare + ^cl . (3.7) 

3.3. Hedgehog meson clouds 

We now examine more carefully the structure of the classical meson cloud. Once 
again a helpful example is the "pions-only" Lagrangian ( |3.1| ). A simpler recasting of the 
Euler-Lagrange equation (|3.2|) comes from nailing the baryon' s center of mass at the origin 
and its (iso)spin orientation at the North pole, X(t) ~ and A{t) = 1, in which case (^3) 
becomes Q 

8V r) 
i-V' + m^X,,(x) + ^ = ^glT.g^S^^) . (3.8) 



from Lorentz invariance (as Jim Hughes has reminded us). So long as the Lagrangian density is 
a Lorentz scalar, meaning that C{x) — > C{x') as x — > x', then L {t) = J d^yiC{x) J d^^C{x') = 
a/i - X^L{t') = (l - + 0{'X^))L{t'), where a/i - X2 is the change-of-frames Jacobian. 
One cannot ignore the Lorentz contraction of the cloud, as the two masses would then be unequal. 

To prove that such terms are suppressed one needs to anticipate the findings of Sec. 3.3 
below, and check that when the classical meson fields are precisely hedgehogs the cross terms 
vanish by symmetry upon spatial integration. They are only nonvanishing to the extent that the 
cloud deviates from the hedgehog ansatz, which it does at a higher order in 1/Nc |^,35|. 



This contrasts with the quantum corrections to the mass of the Skyrmion, which are only 
0{N^). This is because the Skyrmion, unlike the "bare" nucleon, is the solution to an Euler- 
Lagrange equation. 

The reader might be confused about the present definition of mass renormalization, versus 
the use of a conventional mass counterterm 5M — Mbare — Mren. In fact, they are the same thing. 
By definition, 5M must be tuned to cancel the shift in the pole of the propagator away from Mren 
induced by the interactions with the mesons. Since the nonrelativistic propagator i/{kQ — M + ie) 
Fourier transforms to 9{/S.t)e^^^^^^ , this simply means that —i5M/S.t must cancel the meson cloud 
contribution — zMciAt to the effective action. But this condition 5M = — Md is just a rewrite of 
Eq. (3.7), QED. 
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The solution is termed the "static pion cloud" 7fstat(x), whereupon yfci is approximated by 



<l(x,t) = I^iin^(^'))^stat(x) (3.9) 

up to 1/Nc corrections, the primes denoting the baryon's center-of-mass frame as always. 
So it will suffice to focus on Eq. ( |3.8D rather than the more complicated time-dependent 
equation (p.2|). In terms of Tfstat, the quantities Md and Td are explicitly given by 



cl 



j rf^x [i (a.TTstat)' + + y(7rstat) - 3^7^^^^7r,tat • V53(x) ] (3.10) 



and 

Je, = § / d'^T^l,, . (3.11) 



3 

What does Tfgtat look like? An important hint is that the index 'a' lives in zsospace 
on the left-hand side and ordinary space on the right-hand side of Eq. ( p.8|) . Therefore, as 
in the Skyrme model, the solution can generically be found in the maximally symmetric 
"hedgehog ansatz," which equates these two spaces 0: 

^stat(x) = yxF(r) (3.12) 

where r = |x| and x = x/r. In turn, the cloud profile F{r) solves a model-dependent 
nonlinear ordinary differential equation in the radial coordinate r, obtained by plugging 
( CT) into dpD. 

Even if additional mesons are incorporated as per Eq. (|2.11J ), the coupled static equa- 
tions (one for each meson species) are still solved in the hedgehog ansatz, suitably gener- 
alized in the manner familiar from vector-meson-augmented Skyrme models []3T| , ^ : 



Pstat(x) = Ae.afcx'=G(r) , p3f,,(x)=0, astat(x) = Aif(r) , 

(3.13) 



UJ. 



stat 



;x)=0, ^°,,(x)=^i^(r), etc.. 



where 'z' and 'a' label spin and isospin, respectively. Since ~ has been factored 
out explicitly, the profiles F, G, H and K each scale like A*"^. In general they obey model- 
dependent coupled nonlinear radial ODE's. 
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3.4- Ultraviolet divergences and the large-N^ renormalization group 

The above discussion contains something of a cheat: as written, Eq. ( p.8|) only admits 
a mathematicaUy weU-defined solution in the free-pion case V{n) = 0, in which case 
obviously 

^stat(x) oc . (3.14) 

r 

For nonvanishing V, the lack of a solution to this nonlinear equation is due to the exact 
pointlike nature of the 5-function source.0 So, too, the cloud parameters and Xd 
introduced in Eq. ( ^?5|) are actually ill-defined, diverging in the ultraviolet. An ultraviolet 
cutoff is required. This should come as no surprise: in most hadron models the need for such 
a cutoff arises as early in the discussion as the first loop correction in the original Feynman 
graphs (e.g.. Figs. 2b-c). But even if one were to concoct a meson model free of this type 
of divergence at the level of the original graphs, a cutoff A ~ A*"^ on the meson momenta 
would still be required for the self-consistency of the subsequent formalism; specifically, 
one needs to ensure that the baryon always stays in the vicinity of its rest frame even after 
an arbitrary number of interactions with the mesons. In fact, if instead one were to permit 
meson momenta of order Nc, in particular above the NN threshold, it is plausible that 
the entire effective Lagrangian approach breaks down (see Ref. Sec. 8.3). 

The simplest fix, with interesting consequences as we shall see, is to smear out the 
5-function in some manner, 

53(x) 5i(x) , (3.15) 

over a length scale A~^. While details of this regulator should not matter for sufficiently 
large A, we will nevertheless insist that ^^(x) be spherically symmetric, so that the hedge- 
hog ansatz remains valid even for finite A, an enormous technical simplification. 

Now suitably regulated, Eq. ( |3.8| ) (returning to the pions-only example) is easily solved 
numerically for F{r). The solution is shown schematically in Fig. 8. For r well inside the 
classical cloud, whose radius is determined by the parameters of V, the behavior of F{r) 
is highly model-dependent. But for r ^ Eq. ( p.8|) linearizes, and one finds 

F{r) MmiLfl^ + !!hL\e-^-- (3.16) 



This arcane technical point about the lack of a well-defined solution is confirmed in the 
explicit examples of Sees. 5-8 below, as the reader can verify by imagining the following exercise. 
Hold g^'^N fixed (instead of g^j^Sn as in the large-Nc renormalization group) while taking the 
ultraviolet cutoff A — > oo, and confirm that in so doing TTstat has either a singular limit, or no 
limit at all. 
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which follows from ( |3.12 ) and ( |3.14| ). Notice the new parameter, g'!j^^^, which measures 



the height of the exponential tail. While its precise numerical value depends sensitively 
on the choices of V{n) and A, its physical interpretation as the renormalized pion-nucleon 
pseudovector coupling constant is pleasingly model-independent. To our knowledge this 
identification was first made in Sec. 4 of Adkins, Nappi and Witten 0, and has recently 
been confirmed in Ref. with a careful analysis of the LSZ amputation procedure.lll The 
behavior ( p.l6|) is equally valid in the chiral limit: 



o ren 

F{r) -^7^ (massless pions) . (3.17) 

In summary, for any given choice of meson Lagrangian, we have described an explicit 
numerical procedure — Eqs. ( p.7|) and ( p.l6|) — for extracting the physical, renormalized pa- 



rameters Mren, ^ren and Q^^^^pf as fuuctious of the Lagrangian input quantities Mbare, ^^bare, 
g^^% as well as the UV cutoff A. Alternatively, one might wish to fix Mren, ^ren and (7^*^"^, 
say to their experimental values 

Mren + = Mn = 939 MeV , (3.18a) 

O-^ren 
15 

Mren + =Ma = 1232 MeV , (3.186) 

O-^ren 

«™ ^ 2m;7 ■ 

then solve implicitly for Mbare (A), 2^bare(A) and g^^%{A). This latter approach seems the 
most reasonable to us, and will be our philosophy from now on. We call this novel program 
the large-Nc renormalization group, and devote the remainder of this paper to exploring 
its solutions in a variety of illustrative models. 



While both these references discuss this particular issue in the context of the Skyrme model, 
the reader can verify that the conclusions are equally valid for effective Lagrangian models such as 
concern us here, with explicit baryon sources. It should come as no surprise that the LSZ residue 
is sensitive only to the asymptotic behavior ( ^.16| ), as this is a well-known property of Fourier 
transforms. The fact that the renormalized coupling is still pseudovector, like the bare coupling, 
is simply because the Fourier transform of a hedgehog is necessarily proportional to q. The only 
subtle feature is this [0]: the existence of an LSZ pole precisely on the pion mass shell depends 



crucially on a small quadrupole |35] deviation of the Skyrmion (or, in the present context, meson 
cloud) away from the hedgehog ansatz, induced by the baryon's (iso)rotations. This distortion is 
one of the corrections dropped in Eq. (|3.9| ), as it is tangential to our present purposes. 
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4. Free pion Lagrangian and its continuum limit 

Our first example consists simply of free massless pions, 



■^meson — 2 (^M^) ' ('^•■^) 

coupled derivatively to the I = J baryon tower as per Eq. ( |2.8[ ). In the hedgehog ansatz 
( p.l2|) , the static Euler-Lagrange equation ( p.8| ) becomes 



+ 1^' - = 6/-i^7^;^(A)|:5A(r) . (4.2) 

This being a linear equation, it is trivially solved using the method of Green's functions: 

POO o 

F{r) = 6f-'glTA^) dr'r''G{r,r')—SA{r') , (4.3) 

where the massless Green's function that is well behaved at both r = and r = oo is 

G{r,r') = ^ , r<; =min{r, r'j, r> = maxjr, r'} . (4-4) 

3r> 

The renormalized Yukawa coupling g^^^^^ is extracted from the large-distance behavior 
of F as per Eq. ( p.l7|) . With the mild (and relaxable) assumption that 5A(r') has compact 
support, Eqs. ( [4.3|) and imply 



Gg'^TAA) r ^ ^-^f^. (4.5) 

jTrf Jo ^T^J-Kf 

We have made use of an integration by part, plus the requirement that the volume of 
be normalized to unity (regardless of other details of this smearing function). Comparing 
Eqs. ( |i3D and ( pT7| ), we deduce 

glTM = 9':^. (4.6) 

for all A, admittedly not a surprising result for free field theory, but a reassuring sanity 
check on our formalism. 

This lack of any flow is consistent with what might be termed an "exact" Cheshire cat 



picture [|15i[13| (recalling the discussion in Sec. 1). In truth, this is the only model we have 



found where this perfect equality holds. For example, the mildest conceivable modiflcation 
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to the Lagrangian (|4.1|) is to add a pion mass term. In that case the Green's function is, 
instead, 



1 



2mi 



1 \ _ ^ . / 1 . m. 



2 / V 2 



G{r,r') = — ^ ^ _ :::ZL e"^-'^< - + ^ e-"^-'^< + ^ e-"^-'^> (4.7) 



1 m,, 



2 



which properly reduces to ( [4.4| ) as ttItt — * 0, and (|4.6|) is amended shghtly to 



In either variation, massless or massive, the continuum limit A — >^ oo can be safely 
taken, and the "ultraviolet fixed point" that emerges is just what one started with: a 
theory of free pions derivatively coupled to the baryon tower. This is entirely expected, 
since as mentioned earlier V(7f) = is the only case in which Eq. (|3.8| ) as written already 
has a bona fide solution, namely, Eq. ( |3.14| ), and there is no actual need to smear out the 
5-function source. We now turn to more interesting examples where these statements no 
longer hold, and where the breakdown of the "exact" Cheshire cat picture is much more 
severe than Eq. (E 



5. The nonlinear a model and its lack of continuum limit 

For our second example, consider the nonlinear a model for pions, 

Cn^eson = ^Tid^U^d'^U, U = exp {2in ■ T / U) , (5.1) 

again augmented by the bare Yukawa coupling (|2.8[) . The static Euler-Lagrange equation 
( ^78|) now works out to 

F" + - ^sin2F = Qf-'g'^TM^-^S^ir) . (5.2) 
Solving this nonlinear equation for F{r) requires that we specify a smearing of the source. 



For convenience, we follow Ref. and choose a radial step- function 



3Al • r < A-1 
SAir) = i (5-3) 

which is properly normalized to unit volume. The technical advantage, which we exploit 
presently, is that the right-hand side of Eq. ( |5.2|) is now proportional to a true 5-function, 
since 

gS^ir) = -2—S{r-A-') . (5.4) 
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There is also a conceptual advantage: the right-hand side of ( ^.41) means that the baryon 
and meson degrees of freedom only interact at the "bag radius" A~^, which sharpens the 



analogy to the traditional chiral bag |1T0|-[T3 



With this convenient choice of regulator, the prescription for satisfying Eq. (|5.2|) is 
transparent: First solve the homogeneous version of Eq. (|5.2|) , namely 



F" + -F' - ^sin2F = , (5.5) 



for r < ("region I") and for r > ("region II") subject to the boundary condition 
( ^.17| ); next, match the solutions in these two regions, Fi{r) and -Fii(r), at the point 
r = A~^; and finally, read off g^^%{A) from the slope discontinuity, 

^^;^(A) = inUA-'{F;{A-')-F;,{A-')) . (5.6) 

This 3-step graphical procedure is illustrated in Fig. 9. The curves -F„(r) and Fi{r) 
are displayed in Figs. 9a and 9b, respectively. The curve -F„(r) is uniquely specified by 
Qn^NN through the asymptotic formula ( p.l7|) . In contrast, Fi(r) actually stands for an 
entire one-parameter family of curves related by dilatations, -Fi(r) — > Fi(Ar), thanks to the 
scale invariance of Eq. ( ^.5|) . For any specified value of the cutoff A, the scale parameter A 
needs to be tuned so that Fi(A~^) = F„(A~^). Note that Fi, unlike F„, attains a maximum 
value F™^^ = .587r before dropping back down towards n/2. 

Let us discuss the qualitative behavior of the large- A^c renormalization group as the 
cutoff increases from zero. In the infrared regime A <^ if-n-/ Q^nnV^"^ ^ Ao'^ of dir^Ni-^) 
necessarily approaches the free massless pion case, Eq. ( [4.6| ). This is simply because 
the patched-together profile F(r) is small everywhere, in which case Eq. ( |5.2| ) reduces 
to Eq. ( [4.2|) up to 0{F^) corrections. This is precisely the regime studied recently by 
Manohar [jl^, who correctly reproduced the 0{m^) non-analytic correction to the baryon 
mass familiar from one-loop chiral perturbation theory. 

But the behavior of the renormalization group for higher A quickly diverges from 
the free-pions example. Notice that as A passes a first critical value Ai (the point where 
F„(A^^) = 7r/2) a second disconnected solution to Eq. ( ^T2|) emerges, in which F„ intersects 
Fj, not in the branch of the curve labeled "Region A" but rather in "Region B." The flow 
in this branch is still dictated by Eq. ( p.6| ), but since for any given A the value of F/(A~^) 
differs between the two branches, the solutions are distinct. As A increases further, these 
two branches of g^^%{A) gradually approach one another, until at a new critical value A2, 
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defined by -Fii(A2 ^) = F^^^, they coalesce. This latter scale is indicated by a cross in Fig. 9. 
From Fig. 9a one reads off A^^ = .75(3£?;'^;^^/27r/^) ^2 and A2 ^ = .68{3gl^^^/2TTUy/^, or 
in conventional units, 

Ai = 310 MeV, A2 = 340 MeV (5.7) 

using the physical values ( |3.18| ). We calculate 5f^^''^(A2) = ASgl^'^j^. Crucially, for A > A2, 
if one insists that -Fi(O) = 0, then there is no way to match Fi with Fu, hence no solution 
to Eq. 

Thus we have exhibited two "phases" of the model (for want of a better term), the 
first being defined for A G [0,A2], and the second only for A G [Ai, A2]. The critical "bag 
radius" A^^ = .6 fm is the UV scale beyond which the "chiral bag" cannot be formed from 
the 2-derivative pion action ( |5.1| ) alone; higher-derivative terms must be added. Evidently 
this breakdown owes nothing to soft-pion arguments as one is accustomed to — but follows 
solely from large- A^c reasoning. 

To be honest, for A > A2, one can patch together a solution — if one allows the cloud to 
have nonzero winding number (c/. Eq. (6.5a) below). For example, there are four "phases" 
of the model with winding number unity, which are constructed as follows. Leave -Fii(r) 
the same as above, but let Fi{r) Fi{r) + tt, exploiting a discrete symmetry of Eq. ( |5.2|) . 
There are two such solutions, again depending on whether the intersection takes place in 
"Region A" or in "Region B." The two remaining solutions are generated, instead, by the 
discrete symmetry -Fi(r) —Fi{r) + n. Likewise there are four phases of the model for any 
positive winding number n, generated byS Fi{r) ±Fi(r) + nn. Each such phase exists 
only for a finite interval in the cutoff A. In particular — unlike the free-pions example — here 
there is no phase within which one can take the continuum limit A 00. In other words, 
this model lacks an ultraviolet fixed point. 

Of course, there is no reason whatsoever that an effective field theory need have a 
continuum limit. It would be perfectly reasonable to fix A at a finite value less than A2 

Technically speaking, the two independent branches of g^^'^{A), which have coalesced at A2, 
leave the real axis and bifurcate into complex conjugate pairs for A > A2. 

We suspect that these higher winding-number clouds may be, like the hedgehog Skyrmion 
in the B = 2 sector |[l|,|37||, unstable to small deformations in the cloud away from the hedgehog 
ansatz, but we have not looked for any such deformations, neither in the present model nor in the 
ones to follow. 
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where the "chiral bag" stiU makes sense, and to calculate, for example, the static properties 
of this hybrid nucleon, a la Adkins, Nappi and Witten [0. Nevertheless, it is instructive 
to pose the question: In a generic theory, if the continuum limit can in fact be taken, what 
type of UV fixed point might one expect? In the free-meson examples of Sec. 4, g^^'^{A) 
runs to a finite, nonzero value as A — oo. But this behavior must be the exception rather 
than the rule. For, a nonzero limiting value of g^^%{A) suggests that a solution to Eq. ( |37^ ) 
exists even for a non-vanishing right-hand side with an exact 5-function source. Except 
for the free-meson case V = 0, we have yet to discover a differential equation where this 
is possible. Instead, it is far more plausible that g^^'^{A) is the coupling constant of an 
"irrelevant operator," and vanishes in the ultraviolet (as suggested by the factor of A~"^ on 
the right-hand side of Eq. ( |5.6|) ). The resulting continuum theory would then be completely 
independent of the baryonic degrees of freedom, an interesting example of "universality" 
0]. In this event, Eq. ( |3.8| ) admits a solution if (and only if?) the meson Lagrangian 
supports a nontrivial configuration in the absence of a baryonic source — meaning a soliton 
or Skyrmion, either topological or energetic. 

Viewed in this light, it is no surprise that the nonlinear a model coupled to the baryon 
tower has no UV limit, as we have just learned. Plausibly, this is because Eq. ( |5.1|) does 
not by itself support a Skyrmion. The reason is Derrick's famous "no go" theorem (a 
variant of the virial theorem) |T^: if one posits a static Euler-Lagrange solution C/stat(x), 
then the energy functional 

E[U,t.t] = %j d''xTrd,U},,AUst.t , (5.8) 

rather than being stationary, can actually be lowered arbitrarily by a homogeneous rescal- 
ing 

t^stat(x) ^ C/stat(Ax) , E[Ust.t] ^[C^stat] ; (5.9) 

therefore, no such solution can exist.0 There are several known ways to modify the 
nonlinear a model to prevent such a "Derrick collapse." One way, which we examine in 
the following Section, is simply to reduce the dimensionality of space from D = 3 to the 
"critical dimension" D = 2, as Eq. (|5.9|) suggests. As a bonus, the resulting toy model 
turns out to be analytically soluble. Alternatively, and more physically, in Sees. 7 and 8 
we will augment the Lagrangian ( |5.1|) by the 4-derivative "Skyrme term," and explore the 
interesting, and unexpected, consequences. 

This point was apparently missed by Gervais and Sakita [|l4|, whose Eqs. (5.19)-(5.20) admit 
no solution. 
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6. An exactly soluble 2+1 dimensional model with a nontrivial UV fixed point 

Motivated by the above discussion, we would like to construct a 2+1 dimensional 
model which parallels as closely as possible the nonlinear cr model of the previous Section. 
To make the analogy as plain as possible, it is helpful to recall an alternate parametrization 
of the pion field to that given in Eq. ( |5.1|) . Rather than U = exp (2z7f ■ f/f-^), take 

U = uo + iu-T, wo + u^ = l, 7f=-^u, (6.1) 
in terms of which the nonlinear a model Lagrangian (|5.1|) is simply 

^meson = f [(^M^O)' + {^^U^f + {d^U2f + [d^Usf] . (6.2a) 

The natural 2+1 dimensional analog is thenil 

^meson = y [(^M^o)' + {^^n^f + (8^712)'] , nl + nl+nl = l. (6.26) 

Just as the choice of vacuum {uq,u) = (1,0) spontaneously breaks the 0(4) chiral sym- 
metry of (|6.2a|) down to isospin 0(3), so too the vacuum choice (no, n) = (1, 0, 0) in ( |6.2£| ) 



breaks "chiral" 0(3) down to "isospin" 0(2). 
The hedgehog ansatz for the pion cloud, 

(wo,u) = (cosF(r) , xsinF(r)) , (6.3a) 
will prove equally applicable to the lower- dimensional model, 

(no,n) = ( cos F(r) , xsin F(r)) . (6.36) 
Evaluated on these ansatze, the Hamiltonians are quite similar: 

Hsu[F] = I |rf3x(F'^ + 2^) (6.4a) 



versus 



The mesonic sector of this toy model is a simplified version of the "baby Skyrme model" of 
Ref. m. 
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respectively. But there is a major difference in the winding number formulae between the 
two models:lll 

Wsb[F] = - [ d'^^^^^p^ = -[F{0) - F{oo)] (6.5a) 



versus 



/F sin F 
d^x = -i[cosF(0) -cosF(oo)] . (6.56) 

In other words, whereas in three spatial dimensions the hedgehog ansatz is broad enough 
to encompass all integer winding numbers, the same is not true in two dimensions where 
winding number is restricted to the three values { — 1,0, +1} as follows from Eq. (|6.5&|) .B 
This caveat is irrelevant for present purposes, and so we press on. 

There is a well-known rewrite of the 0(3) model in terms of the "conformal variables" 

3| 

ni+in2 _ ni-in2 , „ 

w — , w = . (6.6) 

1 + no 1 + no 

These variables (canonically rescaled by a factor of //2) are closer in spirit to the "old" 

pion representation (|5.1| ), in that there are no extraneous fields such as uq or no that 

need to be eliminated with a spherical constraint. Paralleling the previous Section, we will 

therefore use the w's, not the n's, in constructing the Yukawa coupling. 

Next we turn to the baryons. The obvious "toy" analog of the I = J tower is an 

infinite sequence of states \ v) that transform as (z^, z/), = ±|, ±|, ±|, . . . , under "isospin" 

0{2) = U{1) and spatial U{1) rotations in the x-y plane. A spin x isospin invariant Yukawa 

coupling in the baryon rest frame then has the form 

1^ E ^7. 1^ + 1) (H + H.c. (6.7) 

Here the Qi, are arbitrary complex constants, z = x + iy, and the operator + 1) (z^ | 
simply means that the difference of the U{1) (iso)spin charges between the initial and final 
baryons must be unity. 



23 "Winding number is given, respectively, hy W = J d^:x.eijkTvTiTjTk in three dimen- 

sions, where T^ = U^diU, and W = J :x.tufi n ■ d^n x d^n in two dimensions, where n means 
(no,ni,n2). 

A mild extension of the hedgehog ansatz ( |6.36| ) covers the sectors with winding number n: 
rewrite x as (cos 0, sin 0) and replace this by {cosn6,sinn6) instead. 
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An intelligent way to generate such a coupling is to work in the collective coordinate 
basis \6) ,6 E U{1), analogous to \ A) in 3+1 dimensions. The baryon Lagrangian analogous 
to O is 

'^baryon = " -^bare + |-^bareX^ + l^^bare^^ , (6-8) 

the last term representing free motion on the U{1) manifold. Another of Schulman's path 
integral identities P^ ], 

Kt2) = t 



/ Ve{t) expi dt{- Mbare + ^Xbare^^) = 



(6.9) 



,1 ,3 



with 



Mb^,,, = M^,,e + ^^ and {u\e) = e*^^ (6.10) 

'^-^bare 

equates such motion with the propagation of an infinite tower of energy eigenstates |z/), 
just like Eqs. (PT^)-(PTBD for SU{2). The Yukawa coupling analogous to ( PT^ ) is then 

i^ba.e/^e-^^"(*')52(x') + H.C., (6.11) 

the primed space-time variables referring to the center-of-mass frame of the moving baryon 
as in Sec. 2. Note that Eq. ( |6.11| ) is a special case of Eq. (|6.7D , with all the g^'s equated to 
a single underlying Yukawa constant (/bare (to see this, copy the steps in Eq. ( pj.lO| )). This 
feature too is just like the higher-dimensional example: recall the proportionality rule (v) 
reviewed in Sec. 2.1, and embodied in the matrix elements ( p.lO|) . 

We have assembled all the ingredients necessary to write down the classical Euler- 
Lagrange equation for the meson cloud. The solution, Wc\, automatically sums up all 
contributions to one-meson absorption or emission from the baryon source, to leading order 
in the semiclassical expansion. In solving for Wd, it is convenient as always to boost from 
the Lab frame to the body-fixed frame of the translating, (iso) rotating, Lorentz-contracting 
baryon: 

«;ei(x,t) = e^^(*')«;stat(x') . (6.12) 

In turn, the static meson cloud tt;stat(x) finds a solution in the hedgehog ansatz, composed 
from (|6.3£|) and ( |6?^) : 

z F(r) 

Wstat(x) = - tan——, r = Vii • (6.13) 
r 2 
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Finally, the profile F{r) obeys 
-F'- At sin 2. 

r 2r^ / 2 or 



w/ , 1 w 1 -OF ^bare(A) 2 ^(^) / ^ 

-T H — F r-sin2F = sec -^(>\{r) 



2. (6-14) 



A^<7bare(A) 2 ^(A"^) , 

sec o(r — A 



tt/ 2 

This follows straightforwardly from (|6.4feD , ( |6.11|) and ( |6.13|) , plus an obvious transcription 
of the 5-function regulator ( p.3| ) to two dimensions. 

After this long build-up, we remind the reader of the original motivation behind this 
toy model: by recasting the nonlinear a model in two spatial dimensions, we have side- 
stepped Derrick's theorem, and increased the likelihood of finding a nontrivial UV fixed 
point to the large- A'^c renormalization group equations. We will verify this presently. But 



already our efforts have yielded a bonus: unlike Eq. (5^), the homogeneous variant of 
Eq. ( 6.14 ), namely 

F" + -F' - ^sm2F = , (6.15) 
r 2r^ 

can be solved analytically. Indeed, switching independent variables to logr transforms this 
into the sine-Gordon equation, the solutions to which are the "baby (anti)Skyrmions"il 

F{r) = 2Tan"Vi^ or 7r-2Tan"Vii^ (6.16) 

where fii and /x„ are arbitrary scale constants. 

Figure 10 displays the patched-together solution to Eq. (|6.14|) . In region II (r > 
A~^) we choose the solution that decays to zero, -Fii(r) = n — 2Tan~^//iir, and fix /Xn by 
normalizing the large-r falloff to 5fren/(27r/r) : 

Hu = . (6.17) 

fi'ren 

As before, g^-en is the renormalized Yukawa constant of the model.il Next, pick Fi{r) = 
2Tan~^//ir in region I (r < A~^). Matching Fj to F„ at r = A~^ implies 

A^o 

yren ^ qn 

^: = ^ . (6.18) 



These (anti-) solitons in 2+1 dimensions are also the well-known 0(3) (anti-) instantons in 
two EucUdean dimensions [39|. 

The 'Wc" scalings are gren / ~ V-^c in analogy with the 3+1 dimensional case. Although 
in this model Nc is no longer identified with the number of colors, it still usefully parametrizes 
the semiclassical expansion. To get the factors of n etc. right in the definition of Qren, it suffices 
to solve with a Green's function (cf. Sec. 4) the hnearized version of ( |6T4|) , F" + F'/r - F/r^ = 
iSren / f)S'j^{r), appropriate to the weak-field regime F <^ 1. Alternatively, one can extract the LSZ 
residue of the meson-mass-shell pole as per Ref. H. 
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The running of Qhare is instantly read off from the slope discontinuity in Eq. (|6.14 ) : 



^7ba.e(A) = ^ cos^ (f;(a-1) - f;,{a-': 



^ 2 • (6-19) 



As anticipated, ghaiei-^) approaches Qren in the infrared, while vanishing rapidly (~ A ^) 
in the ultraviolet. 

Note that for any finite value of the cutoff, baryon number in this model is measured 
in the mundane way, by the fermion number operator. Furthermore, the winding number 
( |6.56|) of the patched-together meson cloud is zero, since F{0) = F{oc) = 0. But at infinite 
A the picture looks very different: the baryons have entirely decoupled from the mesons 
thanks to ( |b.l9| ), while the envelope of the sequence of meson clouds with increasing A is 
obviously a soliton or Skyrmion with winding number unity, namely F{r) = Fn^r). In this 
sense, baryon number can be said to have transmuted to winding number at infinite A. 
And in sharp contrast to the previous Section, the UV fixed point of our toy "chiral bag" 
model exists and is nontrivial: a "baby Skyrmion" model. 

As a more sophisticated alternative, one may choose to define baryon number density 
as the sum of the usual explicit fermion number density and the winding number density. 



which is more in the spirit of Refs. |11"12]- For small A, winding number density is 



negligible everywhere; the mundane definition is recaptured, and is entirely concentrated 
inside the bag, r < A~^. But for large A the situation is reversed: explicit fermion number 
density is screened by negative winding number density inside the bag, and the bulk of the 
baryon number is carried, in the form of winding number, by the meson cloud outside the 
bag. In the strict continuum limit, the bag is gone altogether, and only winding number 
remains. This is a greatly simplified variant (with nonrelativistic nucleons, rather than 
spectrally-flowing valence and sea quarks) of the scenario put forward by Goldstone and 
Jaffe 0. 

We now exit toyland, return to 3+1 dimensions, augment the nonlinear a model ( |5.1j ) 
by the well-known "Skyrme term" to overcome Derrick's theorem, and examine under what 
circumstances the statements of the preceding paragraph can, or cannot, be made. 
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7. The Skyrme Lagrangian, with the Adkins-Nappi-Witten value of g. 



ren 
■KNN 



Finally, and most physically, we take for the pion piece of the action the massless 
Skyrme Lagrangian 

>C„ = =^Trd^U^d^U + -^Tr[U^d^U,U^d,U]\ (7.1) 

-LD OZCg 

We continue to assume that the pions are coupled to the I = J tower of explicit baryon 
fields through Eq. (p.8|). Since the two terms in ( |7.1| ) scale oppositely under dilatations 
(P^), Derrick's theorem is avoided, and /^skyrme supports a soliton: the original hedgehog 
Skyrmion (Fig. 11). In their Skyrme-model treatment, Adkins, Nappi and Witten take 
/tt = 129 MeV (vs. 186 MeV experimentally) and Cs = 5.45 in order to fit the nucleon 
and A masses 0. In the present non-Skyrme-model approach, with explicit nucleons and 
Feynman diagrams rather than topology, our above-stated philosophy suggests instead 
that we peg these parameters to their experimental values (although Cg is not very well 
determined by tttt scattering). However, since our present aims are formal rather than 
phenomenological, it is actually best to leave them unspecified. What we do care about 
are their A*";, assignments: 

- ^ - N,. (7.2) 

Thus Nc/h factors out of the action, which justifies our usual semiclassical manipulations. 

Given that the Lagrangian (|7.1j), like the toy model ( |6.2£| ), supports a soliton, what 
might we guess about the large- A'^c renormalization group? Reasoning by analogy with 
the preceding Section, we might expect g^^'^{A) to vanish in the continuum limit, with 
the Skyrmion emerging as the UV fixed point of the family of meson-baryon "chiral bag" 
models. But this naive scenario cannot generally be right! 

To see why not, let us return to a discussion from Sec. 1.4, and think about what 
does — and does not — fiow in our program. What fiows are the bare baryon mass and 
hyperfine mass splitting parameters Mbare(A) and Xbare(A), as well as the bare Yukawa 
couplings g^^%{A), g^f^^{A), etc. What does not fiow are the purely mesonic Lagrangian 
couplings (/tt and Cg in Eq. (|7.1| ), / in Eq. ( |6.2feD , and so forth; see Fig. 12 for a discussion). 



as well as the renormalized parameters Mi-en, ^rem S^nn^ ^^c. These latter quantities, then, 
are independent variables at our disposal in the effective Lagrangian approach — and remain 
so even at the supposed endpoint of the large- A^c renormalization group fiow. Now contrast 
this to the Skyrmion approach. Assume that /tt and Cg have been specified once and for 
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all, and look again at the resulting Skyrmion, Fig. 11. Notice that g'j^^pf, rather than being 
an additional tuneable parameter as we have just argued, is instead fixed by Eq. (|3.17| ) in 
terms of these meson parameters, as originally shown by Adkins, Nappi and Witten 0: 

ISO 

ren ^ anw ~ -^""^ Q^ 

63 Jtt 

The puzzle can now be stated very clearly: If the Skyrme model is in fact connected to 
an effective Lagrangian model by RG flow, then where, why and how has the supposedly 
tuneable Yukawa degree of freedom in the latter disappeared?!! 

The complete resolution of this paradox is the topic of this Section and the next, and 
goes as follows. When g^^^ is tuned to g^%^ precisely, then the Skyrme model does indeed 
emerge as the UV flxed point of the "chiral bag" models, just as in Sec. 6. But for all 
other choices of g^^^^ the physics is closer to that of Sec. 5: the large- A^c renormalization 
group only makes sense up to a critical value of the cutoff Acrit, and admits no continuum 
limit — neither the Skyrme model nor anything else. 

We now flesh this out explicitly. In the hedgehog ansatz (|3.12| ), the Euler-Lagrange 
equation ( p.8|) implied by (|7.1J ) reads 

8 sin^ F\^,, 2 ^, sin 2F / 4 sin^ F 4F'^ 
1 + -^TT^ )F" + -F' — 1 + 



(7.4) 

QA3„bare/'A\ 



We have again adopted the regulator (|5.3| ). As above, we construct solutions Fj{r) and 
F„(r) to the homogeneous variant 

match them up at r = A~^, and extract the running of g^^% from the slope discontinuity: 
Si-W = 1^ ( 1 + ) (f,'(A-') - F,',(A-')) . (7.6) 



ren 

TTNN 



For the remainder of this Section, we focus on the "measure zero" case when g. 
is pegged to its Adkins-Nappi- Witten value ( |7.3D . In that event Fii(r) = -Fskyrmel^), the 



This paradox never came up in the toy model of Sec. 6 because of scale invariance: the 
Skyrmion tt — 2 Tan^^ (47r/r/(/ren) exists for any independently-chosen values of / and Qren- 
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Skyrmion profile of Fig. 11. For -Fi(r), we take tlie family of solutions to Eq. (73) that 
start at the origin and have increasing slope as A itself is increased (Fig. 13). 



We can now evaluate Eq. (|7.6| ). For small r, the Skyrmion profile may be Taylor- 
expanded: 

i^skyrme(?') = TT - CiCgUr + 0{esUr)^ ■ (7.7) 

Therefore, the first term in parentheses in (|7.6| ) is 1 + 8cf + (9(A~^) in the ultraviolet, 
where the numerical slope ci may be read off Fig. 11. The second term in parentheses is 
dominated for large A by 

f;{A-') - ^ , (7.8) 
where C2 is another numerical constant J^l Combining these expressions then gives 

A — 1- oc -| 

glT.W oc — . (7.9) 

Thus, just as in the toy model, the baryons decouple from the pions in the ultraviolet limit, 
at which point ordinary baryon number Yli s {^^^ + A + • • ■ ) effectively transmutes 
into the winding number of the pion cloud, Eq. (|6.5a|) , in the manner discussed at the end 
of Sec. 6. 



8. The Skyrme Lagrangian, with an incorrect choice of g^^j^ 

Finally, we analyze the "generic" case when g^^,^ differs from its Skyrme-model value 
g^%'%- In what way is the solution to Eqs. (|7.4D -([7^) affected? Consider how, in practice, 
one constructs -F„(r) numerically. One first sets F„ ~ 3g^J^^^/{2Trf.j^r'^) in the asymptotic 
regime r ^ (cs/tt)"^, then integrates inwards towards the origin, and in so doing, en- 
counters a surprise: rather than intercepting the y axis at tt as in the previous case with 
9tt^nn = Qn'nN^ diverging for small r as in the nonlinear a model of Sec. 5, F„ invariably 

This behavior is surprising, since generically along curve I, F({r) ~ A instead, as is easily 
argued. However, very near the point where Fi crosses vr, its slope changes over to a A'^ behavior. 
An interesting way of proving this statement (which we first observed numerically) is with the 
calculus of variations. Specifically, if one performs a small deformation of the cloud analogous 
to that shown in Fig. 15 below, and demands that the first variation vanish (as it must), one 
derives Eq. ( |7.8| ), including an explicit expression for C2 in terms of ci and a definite integral of 
the Skyrme Hamiltonian. 
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hits the y axis at a /ia//-integral multiple (n + |)7r, as shown in Fig. 14a. This behavior 
can be confirmed analytically. Linearizing Eq. (^]^) about such half-integral values forces 



F„(r) — ^ [n + l)Ti + const. X r(^+*^)/2 + const, x r^^'^^^/^ 

(8.1) 

+ 0(r3/2) 

The fact that there are two independent constants here, /ii and /i2, shows that this is 
indeed a generic family of solutions to Eq. (73).ii The bizarre oscillatory behavior of this 



essential singularity is verified in Fig. 14b. 

Of course, the family of Fj curves is still given by Fig. 13, and g^^j^%{K) is again read 
off from the slope discontinuity (|7.6| ) at r = A~^, so one can still patch together a bona 
fide solution to Eq. ( [7.4|) . But, as always, this solution is only physically relevant if it is 
locally stable against small deformations of the meson cloud. 

We have carefully investigated this issue of local stability (albeit only within the 
hedgehog ansatz), and identified one dangerous mode, described in Fig. 15. Since we are 
perturbing about an Euler-Lagrange solution, first variations necessarily vanish. Instead, 
our stability analysis focuses on the sign of the coefficient Q of quadratic variations in 
this mode, a positive (negative) value indicating (in) stability. We have calculated (see 
Appendix A for details): 

A — ^oo 

Q — y const. X (e-3^A)^^^ cos log ^2A-i) + 0{K^) . (8.2) 



Therefore, for fixed g^^^j^ 7^ 9^n^j sufficiently large A, the model exhibits alternating 
"phases" of local stability and instability along a large- A''^ renormalization group trajectory. 

The only other self-consistent solutions at small r start at integer multiples of tt, and are of 
the form mr — fj,r + 0{r^), e.g., the Skyrmion itself. Unlike Eq. (|8.l| ), these are only one-parameter 
families (parametrized by /i), so one never "accidentally" stumbles upon them when numerically 
integrating inwards. Amusingly, the same phenomenon holds true at large r as well | |40[ | (cf. 
Fig. 9b): there are generic two-parameter families asymptoting to (n -|- |)7r, and special one- 
parameter families (again, including the Skyrmion) terminating at mr. So the 3-|-l dimensional 
Skyrmion is "special" both at r = and at r = 00, hence has no free parameters. Of course, 
finiteness of the energy requires integer, not half-integer, boundary conditions at both ends, which 
we enforce on our patched-together "chiral bags." 
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Mathematically, this behavior is traceable to the short-distance essential singularities ( |8.1| ) 
in the solutions to Eq. ( |7.5[ ). Each such phase lasts half a period of the sinusoidal oscillation, 
thus: 

unstable phases : 

A G [Acrit , ^Acrit] , A G [K^Aerit , K^Aerit] , A G [ft^crit , ^^Acrit] , " " " (8.3a) 

stable phases : 

A G [ftAcrit , K^Acrit] , A G [ft^Aerit , At'^Aerit] , A G [ft^Acrit , ^^Acrit] , " " " • (8.36) 

Here k = exp(27r/-\/3) is the constant whose logarithm equals half a period, and Acrit is 
defined as the first appearance of this instability as the large- A'^c renormalization group is 
pushed into the ultraviolet. Acrit is analogous to the Landau pole in quantum electrody- 
namics: it is the scale beyond which one cannot push the cutoff while requiring that the 
theory be stable. Of course, the value of Acrit depends sensitively on g^^^^^; in general, the 
closer the latter is to (7^^^, the greater we expect Acrit to be. In this language, the results 
of the previous Section may be understood as the statement that Acrit 00 (stability 
regained) as g'J"^^ g^'^'^. 

As in the nonlinear a model of Sec. 5, one can also patch together pion clouds with 
nonzero winding number, but each of these is again locally stable only for finite ranges in 
A (and potentially locally wTTstable against small deformations away from the hedgehog 
ansatz). While the nature of the mathematical obstruction is somewhat different, the 
conclusion here is the same as for Sec. 5: the Skyrme Lagrangian, coupled to the baryon 
tower, admits no continuum limit when g'tj^^j^ differs from its canonical Skyrme-model value 
g^^^- In this interesting way the "large- A''^ renormalization group as filter" idea is realized. 

We thank Jim Friar, Jim Hughes, Alex Kovner, Aneesh Manohar and Dick Silbar for 
useful comments, and Jessica Binder for help with the artwork. MPM is indebted to the 
Swansea physics department for hospitality during the time that much of this work was 
carried out. This draft was posted on the occasion of ND's 30th birthday. 



Appendix A. Cloud collapse in the Skyrme example when (7^^"^^ 7^ S^nZ- 

**Derivation of the local instability** 

The primary goal of this Appendix is to derive the oscillatory expression ( p.2| ) for 
the Gaussian coefficient Q of the mode shown in Fig. 15. For this purpose, we need to 
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construct the total energy of the cloud, iftot [F] , keeping only the terms of 0{Nc). The usual 
Skyrme-model Hamiltonian density is the sum of the contributions from the 2-derivative 
and 4-derivative terms, 

"^skyrme ~ ^2-deriv 4-deriv 

'■J/„,2 .sin^Fx 1 /sin^F ^F'^sm^F\ (^-1) 



8 V 



Htot is itself the sum of three distinct parts, the contributions from region I, region II, plus 
the negative-definite energy of the Yukawa interaction with the baryonic source: 

iftot[i^] = H,[F] + H,,[F] + ify,kawa[i^] (A.2) 

where 



/•A /"Oo 

H,[F] = 47r / r^rfr T^skyrme , H,,[F] = An / r^rfr Tl^skyrme 

Jo ^A-1 



(A.3) 



and 

i?yukawa[F] = -|/.A^7:;^^^(A)F(A-1) (A.4) 

The point- like form of the latter is due to our special choice of smearing, Eq. ( ^.31) . The 
defining equation for F, Eq. (|7.4|) , is precisely the Euler-Lagrange equation 

^Htot[F] = 0. (A.5) 

Suppose that F{r) is a patched-together solution to Eq. ( [7.4|) , and consider the effect 
of the small deformation of Fig. 15 on each of the three parts in Eq. ( |A.2|) . As constructed, 
the deformation keeps the matching point at r = A~^ always, but changes the value of F 
at this point from Fi(A~^) to — r^)A~^). With this sign convention, a positive value 

of r] lowers the patched-together cloud so that it is closer to the vacuum F{r) = 0, whereas 
a negative value of 77 raises the curve so that it is closer to the Skyrmion of Fig. 11. We 
will treat r] as an infinitesimal parameter, and will find it fruitful to keep all terms through 
order 77'^. The effect on -ffyukawa is immediate^ 

i^yukawa ^ -|/.A^:^-^(A) [f, - ryA-^F/ + ^ry^A-^F/' - ^ry^^A-^F/" + 0{7j^)\ . (A.6) 



From now on, Fi is short for Fi(A ^) and likewise for its derivatives. Note that this equation 
is a perturbative expansion in 77, not in A^^: generically along the solution curve in region I, 
F/(r) ~ 0(A), F/'(r) ~ O(A^), F/"(r) ~ O(A^), and so forth, so that the terms in brackets are 
each nominally 0{A^). (Two important exceptions to this generic behavior along curve I: when 
Fi{r) is near tt then Fl{r) ~ A'^ rather than A as explained earlier, and also, when Fi{r) is near 
it/2 then Fl[{r) ~ A'^^'^ rather than A^ as implied by Eq. (A. 7).) Another possibly confusing 
point: 5^'^^ (A) does not vary under this deformation; it is just a Lagrangian parameter like any 
other for the purposes of the variational calculus. 
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Next, we examine Hj. By design, Fi(r) is supposed to remain a solution to Eq. ( fTSp 
as we perform the deformation. For large A (and therefore, small r), the important terms 
in this equation are those that derive from varying 'H.A-derivi namely: 



-sin -^i(^) w/. X sin2Fi(r) / 4sin F^{r) _ ^ 



The key observation here is that this (approximate) equation is scale invariant, so that the 
solution in region I which satisfies the new perturbed boundary conditions, F"^^(0) = 
and Fj"^*(A-i) = Fj{{l - r])A-^), is just Fj"<^"(r) ^ Fi((l - r])r) for aU r in this range. A 
straightforward Taylor expansion then yields 

.A- 

H, ^ 4n r^drn^.^erw[F,{{l-r])r)] 
Jo 

= M^-V) / r'rfrH4-deriv[Fi(r)] 
= 47r(l - ^) I ^ r^df H4-deriv [F,{r) ] - ^ (2Ff sin^ F, + A^ sin" F^) 

+ 2A2f; sin 2Fi sin^ F, - 2A^ sin" F,^ 
- ^ (f;'(4 - 8 sin^ F,) + 10F;'f;' sin 2F, 

+ A{F[F[" + F/'^) sin^ Fi + A^F^ (12 sin^ Fi - 16 sin" Fi) 

+ (2A2F;'-8A3F;)sin2Fisin2Fi + 6A"sin"Fi) + 0(77") | 

In the first equality we have changed integration variables to f = (1 — 77)r, and exploited 
the scale-covariance of 7Y4-deriv 

This cumbersome expression simplifies considerably if one discards all terms of order 
A^, keeping only those of order A and A^/^. At the point r = A~^ where Fj = Fn, the 
relevant scalings are 

i sin2Fi ^ (n + i)7r - Fj ~ A'^/^ , sin^ Fi ^ 1 , 

(A.9) 

F/j ~ Ai/2 , f; ~ A , F/' ~ A3/2 , F[" ~ A3 
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1 77^ A 

+ 



thanks to Eqs. ( |8.1|) and (|A.7| ). Equation (|A.8|) collapses to 

H, = 47r(l-r7)|^ r'rfr [Fi(r) ] - (2Ff + A^) 

2! 2e2 (^^'^ ^'"^ ^ ^ ^'"^ " ^^^) 

3 A -3 (^-10) 
_ 1. lA^ ( _ 4f;4 + 4F^F^" _ 4a2^/2 _ g^3^/ + 6A^ 

3! 2eg V 

+ 0(77^) + 0{A^) 

Finally we consider region II, for which the deformation at short distances consists of 
adding the constant — r])A~^) — Fi(A~^) to the right-hand side of Eq. (|8.1|) . The 

dominant contribution to the energy in this region comes from the sin'^ -Pii(r) / (2egr^) piece 
of Eq. ( |A.1| ), which is the only term whose integral diverges as A — oo. Another short 
exercise in Taylor expansion givesS 

Hu ^inj^^ __c^^sin4[F„(r) + F:((l-r7)A-i)-F:(A-i)] 

A - 4r7A-iFX + 2r72A-iFf (A.ll) 



e2 



3 A-3 



■ {12 af; f;' + 4f;,f;" -4f;x^) + o{r]^) + c(a°) . 



3! 

The details of the upper limit of integration -Rmax are unimportant, as they do not con- 
ceivably involve divergent terms in A. 

Assembling Eqs. (|A.6| ), ( [A.lOp and ( |A.11| ) yields the series expansion in the small- 



deformations parameter: 

H,ot{v) = Z + TjL + Tj^Q + ri^C + 0{r^^) . (A.12) 

Focusing first on the quantity of greatest interest, the quadratic coefficient Q, we find 

Q = ^(a-2f;'-FF;) sin2F: + 0(AO) , (A.13) 
eg V / 

using Eq. ( |7.6D to eliminate g^'^j^%- Notice that Q is only (9(A^/^), the 0{A) pieces having 
canceled between ( |A.10|) and (|A.11|) . This confirms Eq. ( p.2|) . 



"^^ The appearance of F/j on the right-hand side is due to the observation that J^I^^"^ ^ 



dr \r=A 



r-\ . ^ + 0(A°) precisely when A = (1 ±zV3)/2 as per Eq. 
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The energy of this patched-together cloud is given by the zeroth-order coefficient Z: 



(A.14) 



2n 

= L 



.A-i 

^0 



To ehminate the integral in the second equality, we have used the fact that since the cloud 
is a solution to Eq. (|A.5| ) at ry = 0, the linear term L must vanish: 

»A- 

Numerically, Z turns out to be negative definite for the important case tt, = in Eq. ( |8.1|) , 
corresponding to the coupling range < g^.^^j^ < g^^, as Z is then dominated by the 
attractive Yukawa contribution ( ]A.4D . 
**A nontrivial consistency check** 
Finally, the cubic coefficient is 

471 



C 



3e2A3 

The reason C is interesting is that there is a stringent consistency check on our calculations 
in terms of the ratio R = Q/C, which reads 

F'^ + F''^ + F["^ + = -3A + 0(Ai/2) . (A.17) 
' dF, ' dFI ' dFI' ^ K ) y ) 

Thanks to Eq. (^), this is indeed satisfied by Eqs. ( \KJ^ ) and ( WTWj} . 

To understand the source of this consistency check, look at Fig. 16. It shows that for 
every patched-together solution of the type we have been discussing, there must be another 
nearby solution, labeled "A" and "B", respectively. Algebraically, this follows from the 
cubic polynomial ( [A.12|) , which can be rewritten more explicitly: 

iftot(r7,A) =Z(F,(A-1),f;(A-1),---) + rt(F,(A-i),F;(A-i),---) 

(A.18) 

+ r?3c(F:(A-i),F;(A-i),--.) . 

Aside from the starting solution "A," which is at 77 = by definition, Eq. ( [A.18| ) implies a 
new solution "B" which stationarizes -fftot(^): at ry* = —2R/3 ~ A"^/^. Alternatively, one 
could start one's Taylor expansion at "B," 

Htotiv, A) = ^(F,((l-ry,)A-i), F;((l-r?,)A-i), ■••) 

+ - 77,)A-i) , - r?,)A-i) , ■ ■ ■ ) (A.19) 

+ r?3c(F:((l-r?,)A-'),F;((l-r?*)A-i), •■■) , 
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and demand that the original solution "A" be recovered at ?7 = Consistency between 
these alternative starting points forces 1 — 77** = (1 — t]*)"^. When A ^ 1 so that 77* and 
?7** ^ 1, this in turn implies ry** = — or in other words, 

R{F,{{1 - 77,)A-i) , - 77,)A-i) , ■ ■ ■ ) - -R{F,{A-') , f;{A-') ,■■■), (A.20) 

whereupon Eq. (|A.17| ) follows from a Taylor expansion in 77* . 
**A numerical example** 

Finally we describe a numerical example which bears out the above picture of 
"twinned" solutions "A" and "B," and is further instructive in answering the question, 
What does the locally unstable cloud "A" collapse into? 

Specifically, we have patched together a solution to Eq. j fT^ ) by fixing g^^^^^^ = .Q5g^%^, 
and integrating Skyrme's equation for Fii{r) inwards from infinity, to a matching point at 
r = A~^ where we fix A = 6.5esfn- The "region I" curve -Fi(r) is then the solution to 
Skyrme's equation with boundary conditions -Fi(O) = and Fi(A~^) = Fii(A~^) = .GStt. 
Prom the slope discontinuity one finds g^^% = 8.2e^^/~^ = A6g^^^, while the total cloud 
energy ( [A .21 ) is — 183/7r/es, which is within a few percent of the large-A approximation 
( |A.14D . Since at the matching point F{A~^) > 7r/2, we expect this to correspond to 
the locally unstable configuration "A" in Fig. 16. Can one then uncover a solution with 
F{A~^) < Tx 12, corresponding to the nearby locally stable (but still globally unstable) 
configuration "B"? Sure enough, for the same Lagrangian parameters A and g^^^^i^K), one 
finds a second solution matched at F(A~^) = .4l7r, whose renormalized coupling turns out 
to be tiny, g^^^j^ = .062(7^^^, and whose total cloud energy is — 205/^/es, a little lower 
than that of "A" as Fig. 16 suggests. 

Figure 16 further suggests the existence of a third solution, labeled "C," which is 
energetically much more favorable than either "A" or "B." We can argue for the necessity 
of such a solution, by noticing that for sufficiently negative 77, eventually the cost in gradient 
energy must overwhelm the gain in the Yukawa interaction; therefore there must be a new 
minimum for some finite, negative r\. One naturally guesses that this new solution looks 
more like the Skyrmion, meaning that the value of F(A~^) is closer to tt than 7r/2, and 
that the asymptotics of the tail is closer to Fig. 11. And indeed, for the same Lagrangian 
parameters, we have found a solution for which F(A~^) = I.IStt, g^^^j^ = 1.08(7^5^^, 
the total bag energy is — 338/7r/es. 

While we believe that, for this fixed choice of A, solution "C" is both a local and global 
minimum of the energy, certainly if it is extrapolated far enough into the ultraviolet with 
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the large- A^c renormalization group, it too eventually develops an instability. Although we 
have not pushed the numerics, it is tempting to conjecture that the cloud then collapses 
to one whose renormalized coupling is closer still to g^%^, and that by continuing in this 
way — alternating RG flow with cloud collapse to a new renormalized coupling, followed by 
renewed RG flow to still higher A, and so forth — one iterates one's way to g^%^ precisely. 
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Figure Captions 
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1. Three types of large- A^c models of the strong interactions, and the relationships 



2. (a) A bare meson-baryon coupling (/bare ~ \f^c-, which we shall refer to generically 
as a "Yukawa coupling" even if it involves derivatives. Henceforth, directed lines are 
baryons, undirected lines are mesons, (b) A simple radiative correction to (a). Since 
the 3-meson vertex ~ l/-\/iVc, this graph too ~ Therefore, it is a leading-order 
contribution to (/ren- (c) A more complicated contribution to (/i-en which is likewise leading- 
order. The general rule is: the leading-order graphs are the ones for which, if one erases the 
baryon line(s), one is left with meson tree(s). (d) An example of a subleading contribution 
to (7ren- The purely mesonic loop costs one power of Nc-, so this graph ~ l/\/Nc. 

3. (a) A "seagull" interaction, in which more than one meson interacts with the 
baryon at the same space-time point. The ellipses allow for yet more meson lines than 
shown meeting at the vertex, (b) A Z-graph, in which the baryon runs backwards in time 
over an interval. These arise in the decomposition of Feynman graphs into time-ordered 
diagrams. Time runs upwards in this diagram, (c) An effective seagull, or "Z-guU," 
implied by (b). For Yukawa couplings that obey the It = Jt rule, (c) is suppressed by 

compared to the bare seagull shown in (a). 

4. The complete set of leading-order corrections of the type shown in Figs. 2a-c. The 
oval blob is understood to contain all tree-level meson branchings (no loops). There is an 
explicit sum over the n! attachments of the blob to the baryon line. 

5. The graphical Born-series solution of Eq. ( |3.2| ). The line terminating in a square is 
our notation for tTci, the oval blob contains all tree- level meson branchings, and the Yukawa 
source y is short for the right-hand side of Eq. ( |3.2| ) . 

6. An equivalent rewrite of Fig. 5, with the baryon explicitly drawn (it is implicit in 
the Yukawa source of Fig. 5). The difference with Fig. 4 lies solely in the time-ordering up 
the baryon line (a difference). 

7. (a) The left-hand side of Eq. (p?5|), formed from Eqs. (|2.8| ) and (pTp, in the 
graphical language of Fig. 5. The third summand stands for the sum of all the vertices in 
the potential V{nc\)- Varying (a) with respect to tTci gives Eq. (|3.2|) . (b) Born-series rewrite 
of (a) using the expansion shown in Fig. 5, with combinatoric factors suppressed, (c) The 
meaning of (b) as baryon self-energy and meson-baryon vertex corrections, interpreted in 
the original Feynman-diagrammatic language of Fig. 2. These corrections, shown to the 
right of the baryon line, need to be inserted in all possible locations up this line, as dictated 




> "Chiral 
Effective 



by Eq. (P). 
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8. Typical hedgehog profile F{r). In the large-r regime, all the complicated dependence 
on the meson potential V{Tf) and on the form of the regulator 5a (r) has been reduced to 
the single parameter g^^,^ which measures the height of the tail. 

9. Construction of the patched-together cloud F{r) = 9{A~^ — r)Fi{r) + 6{r — 
A~^)F„(r) for the model of Sec. 5. (a) The curve -F„(r), which is uniquely specified by the 
asymptotic form, Eq. ( |3.17| ). The x axis is in units of the 0{N^) length {3g^^j^/27if-,^y/'^. 
Fii{r) blows up like 1/r for small r. (b) The curve -Fi(r). The cross indicates the curve's 
maximum, ^ .SStt. (We have also marked by a cross in (a) the point where 
-^11 (^) = F^^^ which defines the critical scale r = A^^ discussed in the text.) In Re- 
gion A its slope is positive and in Region B it is negative; eventually the curve asymptotes 
to it/2. The scale of the x axis is purposefully not displayed, because this curve stands for 
the entire one-parameter family of curves related by dilatations: -Fi(r) Fi(Ar). For any 
given value of the cutoff A, this A is to be adjusted so that Fi(A~^) = F„(A~^). 

10. The patched-together solution to Eq. ( |6.14| ), namely: F{f) = 2Tan"^(A2r)6'(A"^- 
f) + (tt — 2Tan~ "^r)^(r — A"-*^), where the dimensionless variables are r = Anfr/g^e-a and 
A = (7j.gnA/47r/. Shown is the curve for A~^ = 4. 

11. The original Skyrmion 0,0]. The x axis is in 0(A^°) length units (esA)"^- 

12. (a) A generic bare 3-meson vertex, (b) and (c): Radiative corrections to the 
vertex, both of which are down by 1/A^c compared to (a), thanks to selection rule (i) 
in Sec. 2.1. Hence, unlike Yukawa couplings, purely mesonic vertices do not run in the 
large- A^c renormalization group at leading order. It is therefore reasonable to fix the bare 
mesonic parameters to the experimental data at the outset. 

13. Family of curves Fi(r) that solve Skyrme's equation. When g^^^^j^ = g^^^, Fu{r) 
is just the Skyrmion profile shown in Fig. 11. 

14. (a) A typical curve -Fii(r) when < g^J^^j^ < f/^^w- explained in the text, rather 
than approaching n for small r, this curve spirals into n/2. The x axis is in units (cs/tt)"^. 
(b) The short-distance behavior of (a) in logarithmic variables, confirming the oscillatory 
behavior of Eq. ( |8.1| ). The variable plotted on the x axis is log(es/7rr), while the y axis is 
{esM-^^HFu - 7r/2) as suggested by Eq. (13). 

15. A dangerous small deformation of the patched-together pion cloud. For alternating 
regions in A, the cloud energy proves to be quadratically unstable against this deformation, 
and the cloud instantly collapses to another configuration closer to the Skyrmion. The 
deformation is defined as follows. Hold the matching point fixed at r = A~^, but raise the 
value of Fi(A~^) = F„(A~^) infinitesimally, while letting the entire curve Fi{r) relax to a 
new solution of Eq. ( [7.5| ) for r < A~^. For r ^ A~^, raise -Fii(r) by a constant amount 
so that the curves stay matched, then merge this curve smoothly to an falloff at large 
distances (the details of this merging are irrelevant for sufficiently large A). 

16. Total cloud energy Htot as a function of the matching-point value Fi(A~^) = 
F„(A~^). For large A, the existence of nearby solutions "A" and "B" follows from the cubic 
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polynomial, Eq. ( |A.12|) , remembering that L = 0, Q = (9(A^/^), and C = 0{A), and that 
the sign of Q depends on whether F(A~^) is greater or less than 7r/2. As A increases further, 
"A" and "B" exchange relative positions every half period in the sinusoidal oscillations of 
Eq. (|8.1|) . Numerical evidence for the global minimum solution "C," nearer the Skyrmion, 
is presented in the text. 
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